OPTICAL SYSTEM FOR MEASURING METABOLISM 
IN A BODY AND IMAGING METHOD 
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BACKGROUND OF THE INVENTION 

The present invention relates to a living 
body optical measurement system and an imaging method 
in the system and more particularly, to a living body 
5 optical measurement system and an imaging method which 
are adapted to measure in vivo information by using 
light and to image results of measurement. 

Desired in clinical medical treatment is a 
system or a method for measuring the interior of a 

10 living body with ease without adversely affecting the 

living body. Measurement using light is very effective 
to the desirability. The first reason for this is that 
the oxygen metabolic function inside the living body 
corresponds to the concentration of a specified pigment 

15 (hemoglobin, cytochrome aa3, myoglobin or the like) in 
the living body, that is, the concentration of a light 
absorber and the concentration of the specified pigment 
can be determined from an absorption amount for light 
(having wavelengths of from visible rays to near 

20 infrared rays ) . The second reason is that light can be 
handled easily by optical fibers. The third reason is 
that optical measurement does not harm the living body 
when used within the safety standards. 

A system which utilizes the advantages of the 

25 living body measurement based on light to irradiate 
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light having wavelengths of from visible rays to near 
infrared rays on a living body and measure the interior 
of the living body from reflection light at a location 
about 10 to 50mm distant from an irradiation position 
5 is described in, for example, patent disclosures of JP- 
A-63-277038 and JP-A-5-300887 . Also, a system for 
measuring CT images of the oxygen metabolic function 
from light transmitting through a living body having a 
thickness of 100 to 200mm, that is, an optical CT 

10 system is described in, for example, patent disclosures 
of JP-A-60-72542 and JP-A-62-231625 . 

Known as a conventional living body optical 
measurement system is an oximeter for measuring the 
degree of oxygen saturation in the artery (JP-A-55- 

15 24004). The oximeter is a system in which light having 
a plurality of wavelengths is irradiated on a living 
body, the transmitting light intensity or reflection 
light intensity from the living body is measured, and 
spectroscopic characteristics of reduced hemoglobin 

20 (Hb) and hemoglobin oxide (Hb02) and pulsation waves 
are utilized to calculate the degree of oxygen 
saturation in the artery. 

Also known as a method for measuring the 
degree of oxygen saturation in tissues of a living body 

25 (average degree of oxygen saturation in both the artery 
system and the vein system) and the hemodynamic amount 
is a method by Jobsus et al (JP-A-57-115232) . This 
method utilizes spectroscopic characteristics of Hb and 
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Hb02 to measure the degree of oxygen saturation and the 
hemodynamic amount in tissues of a living body. 

To add, in the present specification, 
transmitting light, reflection light and scattering 
5 light are not particularly discriminated from each 

other and the intensity of light which is emitted from 
a light source, interacts with a living body and then 
is detected by a photodetector is called the 
transmitting light intensity. 
10 Problems that the invention is to solve are 

three as below. 



(First problem) 

In order to analyze the function of a living 

15 body, a change in hemodynamic movement due to loading 
is sometimes measured from the difference between the 
hemodynamic movement when a load is applied to the 
living body and the change in hemodynamic movement 
during unloading. 

20 The hemodynamic movement during unloading is 

not always constant but changes with time. As an 
example, a time-variable change in transmitting light 
intensity is shown in Fig. 1 which is obtained when 
light is irradiated on a temporal of a subject who lies 

25 quietly on his or her back and the transmitting light 
intensity is measured at a point 3cm distant from a 
light irradiation position. As illustrated in Fig. 1, 
while fluctuation in the measuring system is only about 



0.3%, the living body transmitting light intensity 
changes irregularly and greatly as a whole while 
exhibiting periodical change components. The fluctu- 
ation in transmitting light intensity is attributable 
5 to a change in the hemodynamic movement in the living 
body. 

Even when the subject keeps quiet, the 
irregular change occurs in the transmitting light 
intensity signal in this manner and consequently, when 

10 the transmitting light intensity upon start of measure- 
ment is treated as a reference value, change in hemo- 
dynamic movement due to load is difficult to separate 
from the measured signal. Further, this causes the 
observer to be prevented from deciding whether a time- 

15 variable change in a measured signal displayed on a 

display unit or a time-variable change in hemodynamic 
movement calculated from the measured signal is due to 
fluctuation owned by the living body or is due to the 
application of a load. Accordingly, in the prior arts, 

20 the transmitting light intensity upon start of measure- 
ment is treated as the reference value and therefore, 
there arises a problem that the subject must be kept to 
be quiet to maintain the reference value and the 
measurement cannot be proceeded with for a long period 

25 of time until the signal becomes stable. 



(Second problem) 

It has hitherto been known to optically 
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measure the cerebral cortex under the skull with an 
optical spot by means of a light generating and 
receiving element and a fiber but measurement of an 
image of a hemodynamic state covered with the 
5 protective tissues such as skin tissues and bone 
tissues, that is, measurement with a plurality of 
measuring points has neither disclosed nor suggested. 
For example, when oxygen metabolism changes locally, it 
has been difficult to detect where the change occurs. 

10 For extraction of light signals, the measur- 

ing time, i.e., the number of integral operations of 
measurement must be increased. As a result, the 
measuring time is prolonged and not only a mental 
burden is imposed on the subject but also the operation 

15 efficiency of the system is degraded. 

The present invention intends to solve the 
prior art problems as above. 



(Third problem) 

With the prior arts, the degree of oxygen 
20 saturation in the artery or the hemodynamic movement in 
the living body tissues can be measured. But, the 
prior arts cannot discriminate an change in hemodynamic 
movement due to an overall change in the living body 
from a change in hemodynamic movement due to a local 
25 change in the living body. 

On the other hand, only the change in hemo- 
dynamic movement due to the local change in the living 
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body is sometimes desired to be detected. 

For example, in the cerebrum of the living 
body, a local portion exists which acts in corre- 
spondence to each function of the living body ( herein - 
5 after referred to as a functional portion) and the 

hemodynamic amount or the degree of oxygen saturation 
at the functional portion of the cerebrum changes 
locally in correspondence to an arbitrary function of 
the living body. At that time, if a change in hemo- 

10 dynamic amount or in degree of oxygen saturation at 

only the arbitrary functional portion can be measured 
locally, then the action of the cerebral functional 
portion can be examined in detail, contributing to a 
great importance from the standpoint of medical 

15 science. 

For example, the signal representative of 
fluctuation in transmitting light intensity in Fig. 1 
is difficult to discriminate because an overall 
hemodynamic movement signal in the living body is 
20 accompanied by fluctuation and even when only the local 
hemodynamic movement changes, a signal indicative of 
the change is buried in the fluctuation. 

SUMMARY OF THE INVENTION 

An object of the present invention is to 
25 solve the above problems and to realize a living body 
optical measurement system which images a state of a 
function of the living body by using simplified 
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detectors through measurement within a short period of 
time and a method for imaging results of measurement by 
using the system. 

Another object of the present invention is to 
5 perform measurement which is difficult to achieve with 
the prior arts and in which a local change in hemo- 
dynamic movement is measured separately and discrimi- 
natively from an overall change in hemodynamic movement 
in the living body. 

10 (For the first problem) 

When a change in hemodynamic movement due to 
load is measured, measurement is carried out by 
alternately providing time during which the load is not 
applied to a living body (unloading time) and time 

15 during which the load is applied to the living body 

(loading time). Here, where a signal measured by the 
living body optical measurement system (measured 
signal) is Sra(t), a signal due to a change in hemo- 
dynamic movement during unloading (non-load signal) is 

20 Str(t), and a signal due to a change in hemodynamic 
movement during loading (load signal) is Sl(t), the 
measured signal Sm(t) is given by the following 
equation ( 1 ) : 

Sm(t) = Str(t) + Sl(t) (1), 

25 t being the measuring time. 
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In the present invention, a signal during the 
non-load time is extracted from the measured signal 
Sm(t) to predict the function Str(t) indicative of the 
non-load signal (estimated non-load signal) and the 
5 load signal Sl(t) is determined from the difference 
between the measured signal Sm(t) and the estimated 
non-load signal Str(t). Further, by displaying the 
determined measured signal and the predicted non-load 
signal at a time, thus making it easy to decide whether 

10 fluctuation in the measured signal is due to fluctu- 
ation in the load or due to fluctuation attributable to 
the living body during unloading. 

Determination of the function Str(t) can be 
effected by inputting an arbitrary function having 

15 indefinite coefficients into a computer through, for 
example, a keyboard and determining the indefinite 
coefficients by, for example, the method of least 
squares such that the function fits optimally to the 
non-load signal. The load signal Sl(t) does not fall 

20 to zero as soon as the load is removed from the living 
body and therefore, by setting the predetermined 
relaxation time following the loading time and deter- 
mining the function Str(t) by using the measuring time 
corresponding to the unloading time exclusive of the 

25 relaxation time, a more accurate function Str(t) can be 
determined . 

The above function Str(t) can be determined 
such that a single function can cover a plurality of 



load times, for example, the entire measuring time or 
can be determined every load time so as to cover only 
the respective load times . With a method of 
determining functions St(t) for the respective load 
5 times by using measured signals Sm(t) obtained before 
and after each load time, high estimated accuracy can 
be obtained. 

(For the second problem) 

To achieve the above object, a living body 

10 optical measurement system according to the present 
invention comprises a plurality of light irradiation 
units for irradiating light rays having the wavelength 
range of from visible rays to near infrared rays on a 
subject, a plurality of light receiving units for 

15 detecting light rays irradiated from the light irradi- 
ation units and transmitting through the interior of 
the subject, a memory for storing, in time sequence, 
signals detected by the light receiving units and 
delivered out of the respective ones of the plurality 

20 of light receiving units, an arithmetic unit for 

performing conversion into signals of measured objects 
at a plurality of measuring points by using the signals 
stored in the memory, and an image preparation unit for 
determining a signal at a measuring position from the 

25 output of the arithmetic unit is measured presumptively 
and displaying the determined signal as an image 
indicative of an intensity signal on a two-dimensional 
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display screen. In particular, each of the plurality 
of light irradiation units includes a plurality of 
light sources having different wavelengths, modulators 
for modulating the light rays of the plurality of light 
5 sources with different frequencies, and wave guides for 
guiding a plurality of modulated light rays to the 
irradiation positions, and each of the plurality of 
light receiving units includes a splitter for splitting 
the intensity of light from each of the plurality of 
10 light sources having the different wavelengths. 

The presumptive measuring position location 
is intermediate between the light irradiation unit and 
the light receiving unit and more specifically, the 
exact presumptive measuring position is an intermediate 
15 portion between a living body surface position irradi- 
ated with the light form the light irradiation unit 
and a light receiving surface position of the living 
body. Here, however, the light irradiation unit is 
close to the light receiving unit and therefore, it is 
20 not particularly problematic that the presumptive 

measuring position can be replaced with substantial 
half the distance between the center of the light 
irradiation unit and the light receiving unit. 
Further, in one method for determination of the inter- 
25 mediate position, light from the light irradiation unit 
is irradiated on the living body and signals are 
detected by two light receiving units disposed at sites 
which are symmetrical to the light irradiation unit. 
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The difference signal between these signals is produced 
and an intermediate location between the light receiv- 
ing unit location and the light irradiation unit which 
is obtained by positioning the light receiving units 
5 through such an adjustment of the two light receiving 
units that the difference signal is rendered to be zero 
level . 

To make up the above description, when the 
function of the interior of a living body is measured 

10 using the above living body optical measurement system, 
light irradiation positions of the plurality of light 
irradiation means are distributively arranged on a 
measuring portion of a subject, a plurality of light 
receiving portions of the plurality of light receiving 

15 means are disposed around the distributively arranged 
light irradiation positions, respectively, and for 
light signals detected by the plurality of light 
receiving means, a position on a perpendicular vertical 
to the living body surface, which position is an 

20 intermediate point between the light irradiation posi- 
tion and the light detection position and extending to 
the interior of the living body is set as a presumptive 
measuring point. The signal at the presumptive 
measuring point is calculated by the light signal 

25 detected by the plural light receiving means. This is 
because according to spatial characteristics of light 
density irradiated from the light irradiation position 
and then reaching the detection position, the density 
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is high immediately below the irradiation position and 
the detection position near the surface by virtue of 
scattering of light from the surface but when an 
arbitrary depth determined by the distance between 
5 light irradiation position and detection position and 
the light scattering characteristics of the living body 
is exceeded, the density becomes the highest even at an 
intermediate location between the irradiation position 
and the detection position and eventually, the sensi- 

10 tivity becomes the highest at the intermediate loca- 
tion. Light signal intensity detected in correspond- 
ence to the aforementioned presumptive measuring point 
and measuring point is displayed on a two-dimensional 
image. When the light signal intensity is desired to 

15 be displayed as a topography image, signals at 

unmeasured locations can be obtained in the form of 
interpolated signals associated with the aforementioned 
presumptive measuring points . 

When the above subject is a living body, the 

20 distance between the irradiation position and the 

detection position is preferably 10 to 50mm. Here, 
determination of the factor for determining the maximum 
distance has relation to the intensity of irradiated 
light and attenuation in the living body. 



25 



(For the third problem) 

Advantageously, in the present invention, 
light rays are irradiated from a desired single site 
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or a plurality of sites on a living body, two sites of 
a detection position at which a local change is 
measured as a change in signal and a detection position 
at which a local change is not measured as a change in 
5 signal are set to be substantially equidistant from a 
light irradiation position, transmitting light inten- 
sity levels are detected at the respective detection 
positions, and the difference between transmitting 
light intensity levels at the two sites is produced, so 

10 that a fluctuation component in the living body common 
to the two detection positions can be removed to permit 
a slight change in one of the light receiving units to 
be detected with high sensitivity. When the afore- 
mentioned two sites of detection position cannot be 

15 found, the light irradiation position is displaced to 
find the detection positions. Namely, like a stetho- 
scope, wanted measuring positions can be searched. 

Preferably, transmitting light rays are 
received at two sites of detection positions 

20 equidistant from the incident position and differently 
positioned, transmitting light intensity levels at the 
respective detection positions are converted into 
electric signals by using photoelectric conversion 
elements such as photodiodes or photomultiplier tubes 

25 (hereinafter, electric signals meaning the transmitting 
light intensities will be referred to as transmitting 
light intensity signals ) , the individual transmitting 
light signal intensity levels are subjected to 
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logarithmic conversion by means of logarithmic 
amplifiers, and a transmitting light intensity signal 
at the first detection position and a transmitting 
light intensity signal at the second detection position 
5 are then amplified and detected by differential 
amplifiers . 

By modulating light emitted from the light 
source in intensity and extracting only a frequency 
component for intensity modulation of a detected 
10 signal, noises due to external disturbance can be 
removed . 

The light source can be connected to the 
light irradiation position through an optical fiber and 
the light detection position can be connected to the 

15 photodetector through an optical fiber. 

Since, in the present invention, information 
about a measuring position can be determined substan- 
tially definitely by the position at which light is 
irradiated on a subject by the light irradiation means 

20 and the position of the light receiving means, the 

signal processing for displaying the information as an 
image can be conducted easily at a high speed. Also, 
with the light receiving means disposed about 10 to 
50mm closely to the light irradiation position, 

25 transmitting light is utilized to obtain the detection 
intensity which is sufficiently high, amounting to 
about 6 order or more higher than light transmitting 
through a living body of about 100 to 200mm. 
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Therefore, measurement can be carried out with simpli- 
fied photodetectors and can be completed within a short 
period of time. 

For example, when the object to be 'measure 
5 (subject) is the head, it is known as reported in, for 
example, "Intracerebral penetration of infrared light" 
by Pa trie W.McCormic el al, Journal of Neurosurgery 
published in Feb., 1992, Vol.76, paragraphs 315-318) 
that for the distance between the irradiation position 

10 and the detection position being at least 30mm, detec- 
tion light transmits through the skin and skull to 
reach the surface portion of the cerebrum, i.e., the 
cerebral cortex. It is also known from characteristics 
of light propagation in the living body that informa- 

15 tion from a location which lies on a perpendicular 
vertical to the surface of the living body and 
extending to the interior of the living body from an 
intermediate point between the irradiation and detec- 
tion positions is contained most richly in light 

20 detected at that detection position. Such character- 
istics are reported in, for example, "Monte Carlo 
simulation of photon path distribution in multiple 
scattering media" by Shecha Feng et al, SPIE, 
Proceedings of photon migration and imaging in random 

25 media and tissues, Vol. 1888, paragraphs 78-89, (1993). 

The living body optical measurement system of 
the present invention needs a number of light irradia- 
tion means and a number of light receiving means for 
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measurement at many positions but as will be described 
in embodiments to be described later, the system is 
effective to measurement at a partial position and 
images can be obtained through a simplified arithmetic 
5 processing in which measurement results obtained at a 
plurality of measuring points are interpolated in 
association with the respective measuring points. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a graph showing fluctuation in 
10 living body transmitting light intensity obtained with 
a conventional system. 

Fig. 2 is a block diagram showing the 
construction of an embodiment of a living body optical 
measurement system according to the present invention. 
!5 Fig. 3 is a diagram for explaining an 

embodiment of an imaging method using the measuring 
system of Fig. 2. 

Fig. 4 is a graph depicting a time-variable 
change in measured signal at a measuring point in the 
20 above embodiment and a time-variable change in 

estimated non-load signal 15 determined from the 
measured signal. 

Fig. 5 is a graph depicting a time-variable 
change in relative change amount of hemoglobin 
25 concentration at a measuring point in the above 
embodiment . 

Fig. 6 is a graph depicting a topography 
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image in the above embodiment. 

Fig. 7 is a graph depicting another example 
of the topography image in the above embodiment. 

Fig. 8 is a graph depicting an example of 
5 display of topography image in the above embodiment. 

Fig. 9 is a diagram for explaining a method 
for coordinate conversion in another embodiment of the 
living body optical measurement system of the present 
invention. 

10 Fig. 10 is a graph showing an example of 

display in the measuring system of the present 
invention . 

Figs. 11A and 11B are graphs showing examples 
of display according to the measuring system of the 
15 present invention. 

Figs. 12A and 12B are graphs showing examples 
of display according to the measuring system of the 
present invention . 

Fig. 13 is a graph showing an example of 
20 display according to the measuring system of the 
present invention. 

Fig. 14 is a graph showing an example of 
display according to the measuring system of the 
present invention. 
25 Fig. 15 is a block diagram for explaining the 

construction of a system according to another embodi- 
ment of the present invention. 

Fig. 16 is a block diagram for explaining the 
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construction of a system according to another embodi- 
ment of the present invention. 

Figs. 17A and 17B are schematic sectional and 
bottom views of a light detection probe in another 
embodiment of the present invention. 

Fig. 18 is a schematic bottom view showing 
another embodiment of the light detection probe. 

Fig. 19 is a schematic diagram for explaining 
an example of use of the light detection probe. 

Fig. 20 is a block for explaining the 
construction of a system according to another embodi- 
ment of the present invention. 

Fig. 21 is a block for explaining the 
construction of a . system according to another embodi- 
ment of the present invention. 

Fig. 22 is block for explaining the construc- 
tion of a system according to another embodiment of the 
present invention. 

Fig. 23 is a block for explaining the 
construction of a system according to another embodi- 
ment of the present invention. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 
(Embodiment 1) 

Embodiments of the present invention will be 
described hereinafter. 

Fig. 2 shows the construction of an embodiment 
of a living body optical measurement system according 
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to the present invention. The present embodiment is an 
example in which the living body optical measurement 
system is applied to the measurement of a change in 
hemodynamic movement attributable to a cerebral 
5 function (a relative change amount of concentration of 
hemoglobin oxide and reduced hemoglobin). A specified 
portion of the cerebrum is related to control of an in 
vivo specified function (for example, moving a part of 
a body such as fingers) and the hemodynamic movement at 

10 the specified cerebral portion is changed by activating 
the specified function. In one form of the present 
embodiment, the living body optical measurement system 
can also be used to measure changes in hemodynamic 
movement under the application of a load for activation 

15 of the aforementioned specified function, for example, 
moving fingers and display measured changes in the form 
of a contour map on a two-dimensional plane image 
indicative of the cerebral portion. 

As shown in Fig. 2, there are provided in the 

20 present embodiment a plurality of light sources 2a to 
2d having different wavelengths (the light sources 2a, 
2c and the light sources 2b, 2d having the same wave- 
lengths, respectively, which are in the range of from 
visible rays to near infrared rays ) , modulators for 

25 modulating the intensity of light rays of the plurality 
of light sources 2a and 2b (2c and 2c) by means of 
oscillators la and lb (lc and Id) having mutually 
different frequencies, a plurality of light irradiation 
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means for irradiating light rays from couplers 4a (4b), 
adapted to couple an intensity-modulated light ray 
propagating through an optical fiber 3a ( 3c ) and that 
propagating through an optical fiber 3b (3d), onto 
5 different positions on the scalp of a subject 6 

standing for an object to be examined through the 
medium of optical fibers 5a (5b), and a plurality of 
light receiving means comprised of photodetectors 8a to 
8f provided for a plurality of optical fibers 7a to 7d 

10 for light detection and a plurality of light detection 
optical fibers 7e and 7f having their tip ends posi- 
tioned near the light irradiation positions of the 
plurality of light irradiation means at locations 
equidistant (assumed herein to be 30mm) from the light 

15 irradiation positions. Living body transmitting light 
rays are collected to optical fibers by means of the 
six light detection optical fibers 7a to 7f and are 
photoelectrically converted by the photodetectors 8a 
to 8f , respectively. The light receiving means is 

20 operable to detect and convert light reflected inside 

the subject into an electric signal and a photoelectric 
conversion element represented by a photomultiplier 
tube or a photodiode is used as the photodetector 8. 

Electric signals indicative of the living 

25 body transmitting light intensity levels which are 

photoelectrically converted by the photodetectors 8a to 
8f (hereinafter referred to as living body transmitting 
light intensity signals) are inputted to lock-in 
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amplifiers 9a to 9h, respectively. Since the photo- 
detectors 8c and 8d detects living body transmitting 
light intensity levels collected by the light detection 
optical fibers 7c and 7d which are equidistant from 
5 both of the optical fibers 5a and 5b, signals from the 
photodetectors 8c and 8d are split into two systems so 
as to be inputted to the lock-in amplifiers 9c and 9e 
and the lock-in amplifiers 9d and 9f . The intensity 
modulation frequencies from the oscillators la and lb 

10 are inputted, as reference frequencies, to the lock-in 
amplifiers 9a to 9d and the intensity modulation 
frequencies from the oscillators lc and Id are 
inputted, as reference frequencies, to the lock-in 
amplifiers 9e to 9h. Accordingly, living body trans- 

15 mitting light intensity signals associated with the 

light sources 2a and 2b are separately delivered out of 
the lock-in amplifiers 9a to 9d and living body trans- 
mitting light intensity signals associated with the 
light sources 2c and 2d are separately delivered out of 

20 the lock-in amplifiers 9e to 9h. 

Exemplarily, for contour map display, the 
separated transmitting light intensity signals of 
individual wavelengths delivered out of the lock-in 
amplifiers 9e to 9h are subjected to analog to digital 

25 conversion by an analog to digital converter 10 and are 
then stored in a memory 12 provided internally or 
externally of a computer 11. During or after measure- 
ment, the computer 11 uses the transmitting light 
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intensity signals stored in the memory to calculate 
relative change amounts of concentration values of 
hemoglobin oxide and reduced hemoglobin which are to be 
determined from detection signals at individual 
5 detection points and stores the calculated amounts in 
the memory 12 as time-variable information at a 
plurality of measuring points m. The above calculation 
will be described later in greater detail. A display 
controller 30 converts the signals stored in the memory 

10 means 12 into display signals for a display unit 13 

such as a CRT and displays them on the display unit 13. 
In the display signals, measuring positions are 
converted into coordinates on the display plane of the 
subject and treated as signals for intensity signal 

15 (relative change amounts of concentration values of 
hemoglobin oxide and reduced hemoglobin) contour map 
display at the coordinate positions. 

By using the living body optical measurement 
system according to the present embodiment, relative 

20 change amounts of concentration values of hemoglobin 

oxide and reduced hemoglobin in the living body can be 
measured easily at a high speed. The expanded con- 
struction in which the number of light incident points 
(light irradiation positions) and the number of light 

25 detection points are increased can be feasible with 

ease by increasing the number of intensity modulation 
frequencies of the light sources, of the light sources, 
of the photodetectors and of the lock-in amplifiers. 
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With the present living body optical measurement system 
used, the spectroscopic positions and the light 
irradiation positions can be separated in accordance 
with the intensity modulation frequencies and there- 
5 fore, even when the number of the light irradiation 

positions is increased, it suffices that the number of 
wavelengths of irradiation light rays at the respective 
light irradiation positions equals the number of 
absorbers to be measured, and the wavelength of an 

10 irradiation light ray need not particularly be changed 
for the respective light irradiation positions. 
Accordingly, the number of wavelengths of irradiation 
light rays used is small and an error due to the 
influence of scattering which changes with the 

15 wavelength can be decreased. 

Figure 3 is a diagram for explaining an 
embodiment of an imaging method according to the 
present invention which uses the living body optical 
measurement system, showing the relation among the 

20 light incident point, the light detection point and the 
measuring point in the above method. The imaging 
method of the present embodiment is a method for 
preparing images of relative change amounts of con- 
centration values of hemoglobin oxide and reduced 

25 hemoglobin at the head of a subject, wherein four 

incident points and four detection points are provided 
at the left temporal participating in the motion 
function of right fingers of the subject so as to 
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measure living body transmitting light intensity levels 
and results of measurement obtained under the applica- 
tion of a load of right-finger motion and a load of 
left-finger motion are imaged. 
5 As shown in Fig. 3, light incident points 17a 

to 17d and detection points 18a to 18d are disposed on 
the left temporal of a subject 16. Here, the respec- 
tive light incident points are in correspondence 
relationship with the respective detection points 

10 through tens sets of 17a-18a, 17a-18b, 17b-18a, 17b- 
18b, 17b-18c, 17b-18d, 17c-18b, 17c-18c, 17d-18c and 
17d-18d. The distance between the corresponding light 
incident point and detection point is 30mm. Further, 
as described in the previously-described "Monte Carlo 

15 simulations of photon path distribution in multiple 

scattering media" by Shechao Feng et al, time-variable 
changes in the relative change amounts of concentration 
values of hemoglobin oxide and reduced hemoglobin to be 
determined from measured signals at the respective 

20 detection points best reflect information from an 

intermediate point between the corresponding incident 
point and detection point and hence, presumptive 
measuring points 19a to 19 j are each set in the middle 
of the correspondence relation between each incident 

25 point and each detection point. Information at the 

presumptive measuring points 19a to 19 j is determined 
and the magnitude of the information is displayed in 
the form of a contour map, a light and shade map or a 
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color discrimination map on the two-dimensional plane 
as shown in Fig. 3. 

Next, there will be described an embodiment 
of a method according to the present invention for 
5 determining relative change amounts of concentration 

values of hemoglobin, that is, changes in concentration 
values of hemoglobin at a specified cerebral portion 
obtained when a specified in vivo function (for 
example, moving a part of body such as fingers) is 

10 activated, from measured signals at the respective 
light detection points. 

Fig. 4 is a graph illustrating a measured 
signal 14 at one of the detection points 18a to 18d of 
the living body optical measurement system and a time- 

15 variable change in estimated non-load signal 15 deter- 
mined from the measured signal 14 in the embodiment of 
Fig. 3. In the graph, abscissa represents measuring 
time and ordinate represents relative concentration 
change amount. 

20 The estimated non-load signal 15 is 

determined by removing from the measured signal 14 
signals occurring during time Tt for a load to be 
applied (loading time) and during time T2 for the 
signal to recover its intact form following loading 

25 (relaxation time) and fitting an arbitrary function to 
the measured signal 14 occurring during load preceding 
time Tl and load succeeding time T3 through the use of 
the method of least squares. The present embodiment is 
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handled by using a quadratic linear polynomial as the 
arbitrary function and setting respective times to 
Tl=40 seconds, T2=30 seconds, Tt=30 seconds and T3=30 
seconds . 

5 Fig. 5 is a graph illustrating time-variable 

changes in relative change amounts of concentration 
values of hemoglobin oxide and reduced hemoglobin 
(hereinafter represented by ACoxy(t) signal 20 and 
ACdeoxy(t) signal 21, respectively) at one measuring 

10 point. In the graph, abscissa represents measuring 
time and ordinate represents relative concentration 
change amounts. A hatched time interval corresponds to 
load applying time (motion period of right fingers). 
From the measured signal 14 and estimated non-load 

15 signal 15 of two wavelengths indicated in Fig. 4, the 

above relative change amounts of concentration value of 
hemoglobin oxide and reduced hemoglobin (Hb02, Hb) 
under the application of a load are determined through 
the following arithmetic operation processing. 

20 For a wavelength A, the relation between 

estimated non-load signal Str(A, t) and light source 
intensity 10(A) is given by the following equation (2) 
by separating in vivo light attenuation into scattering 
and absorption: 

25 -Ln{Str(A, t)/l0(A)} 

=€oxy(X). Coxy(t). d+6deoxy(X). Cdeoxy(t). d+A(A)+S( A.) 

(2) 

where 



€oxy(A): extinction coefficient of hemoglobin 

oxide at wavelength A. 
€deoxy(A): extinction coefficient of reduced 

hemoglobin at wavelength A 
A (A.): attenuation due to absorption by other 

substances than hemoglobin at wavelength A 
S ( A ) : attenuation due to scattering at 

wavelength A 
Coxy(t): concentration of hemoglobin oxide at 

measuring time t 
Cdeoxy: concentration of reduced hemoglobin at 

measuring time t 
d: in vivo effective optical path length (in a 
region of interest) 

Also, the relation between measured signal 
Sm(A r t) and light source intensity 10(A) is given by 
the following equation ( 3 ) : 
-Ln{Sm(A,t)/lO(A)} 

=€oxy(A). {Coxy(t)+C'oxy(t)+Noxy(t)}. d 

+Gdeoxy ( A ) . {Cdeoxy ( t ) +C ' deoxy ( t ) +Ndeoxy ( t ) >. d+A' ( A ) + 
S'(A) ..(3) 
where 

C'oxy(t): change in concentration of hemoglobin 
oxide under the application of load at 
measuring time t 

C'deoxy(t): change in concentration of reduced 
hemoglobin under the application of 
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load at measuring time t 



Noxy(t): noise or high frequency fluctuation in 



concentration of hemoglobin oxide at 



measuring time t 



5 



Ndeoxy(t): noise or high frequency fluctuation in 



concentration of reduced hemoglobin at 



measuring time t 



Here, if A (A,) and S(A) remain unchanged under loading 
and unloading, that is, if a change in the measured 
10 signal under the application of load is due to only 

changes in concentration values of hemoglobin oxide and 
reduced hemoglobin, the difference between equations 
(2) and (3) is given by the following equation (4): 

Ln{Str(X,t)/Sm(A,t)}=eoxy(A.) {C'oxy(t)+Noxy(t)>d 



change amounts of concentration values of hemoglobin 
oxide and reduced hemoglobin under the application of 
load are represented by ACoxy(t) and ACdeoxy(t) and 
20 defined by the following equations 



15 



+ edeoxy ( X ) { C ' deoxy ( t ) +Ndeoxy ( t ) } d 



. . . . (4) 



Here, time-variable changes in relative 



ACoxy(t)={C'oxy(t)+Nox(t)}d 



ACdeoxy ( t ) = { C ' deoxy ( t ) +Ndeoxy ( t ) } d. 



(5) 



Usually, it is difficult to specify d and 
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therefore, the dimension of the concentration change 
amount is herein the product of concentration and 
distance d. 

But, in equation (5), distance d acts equally 
5 on ACoxy and ACdoxy and so equation (5) is considered 
as indicating relative change amounts of hemoglobin 
concentration values . When two wavelengths are used 
for measurement, obtained equation (4) is reduced to a 
simultaneous equation with two unknowns for ACoxy(t) 

10 and ACdeoxy(t), so that ACoxy(t) and ACdeoxy(t) can be 
determined from estimated non-load signal Str(X, t) and 
measured signal Sm(X, t) for each wavelength. Further, 
what is indicated by ACoxy(t) and ACdeoxy(t) during a 
period other than loading time and relaxation time can 

15 be formulated by C'oxy(t)=0 and C'deoxy(t)=0, which 
indicate noises or high frequency fluctuations, 
attributable to the living body, of concentration of 
hemoglobin oxide and reduced hemoglobin. The above 
procedure is carried out for 0 to 140 seconds to obtain 

20 ACoxy(t) signal 20 and ACdeoxy(t) signal 21 of Fig. 5. 

Figures 6 and 7 show contour map images 
(topographic images) prepared from time-variable 
changes in relative change amounts of concentration 
values of hemoglobin oxide at the respective measuring 

25 points under the application of loads of left-finger 
motion and right-finger motion of a subject, respec- 
tively. A method for preparation of the topography 
images is such that a time integral (alternatively. 
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time average) of relative change amount ACoxy(t) signal 
20 during load applying time (hatched period in Fig. 5) 
is calculated by the computer 11 and values between 
individual measuring points are determined through 
5 linear interpolation in X-axis and Y-axis directions. 
In addition to the contour map, the topographic image 
may include a monochromatic light and shade image and a 
discrimination display in color. It will be seen by 
comparing the images of Figs. 6 and 7 that the con- 

10 centration of hemoglobin oxide clearly increases at 

specified positions during the right hand motion. By 
displaying this type of spatial distribution informa- 
tion as images, rapid and easy recognition of measure- 
ment results can be ensured. While the images shown in 

15 Figs. 6 and 7 are prepared from the time integral 

values of concentration relative change amounts during 
the load applying time, a topography image can also be 
prepared similarly from relative change amounts of 
concentration values of hemoglobin oxide which are 

20 measured every constant measuring time at the respec- 
tive measuring points. When a plurality of topographic 
images thus prepared in accordance with order of the 
measuring time points or displayed as a moving picture, 
time-variable changes in the relative change amounts of 

25 concentration values of hemoglobin oxide can be 
obtained . 

Further, by calculating time-variable changes 
in relative change amounts of concentration values of 



hemoglobin oxide at an arbitrary measuring point as 
well as a self correlation function and a mutual 
correlation function of time-variable changes in 
relative change amounts of concentration values of 
5 hemoglobin oxide at that arbitrary measuring point and 
another measuring point, a topographic image can be 
prepared from the correlation functions at the 
individual measuring points. The correlation function 
at each measuring point is a function defined by time 

10 shift x and therefore, by preparing topography from a 
value of the correlation function at the same time 
shift x and displaying the topography in accordance 
with order of x or displaying it in the form of a 
moving picture, hemodynamic movement which changes by 

15 participating in activation of the cerebral function 
can be visualized. 

Here, the relative change amount of concen- 
tration of hemoglobin oxide is typically used for 
explanation but the relative change amount of concen- 

20 tration of reduced hemoglobin or the total hemoglobin 
concentration change amount which is calculated by 
adding the relative change amount of concentration of 
hemoglobin oxide and the relative change amount of 
concentration of reduced hemoglobin may be used to 

25 prepare topographic in a similar way. 

Figure 8 shows a display example in which a 
topography image 22 prepared through the above- 
described method is superimposed on a cerebral surface 
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image 23. Since the topographic image 22 is illus- 
trative of a change in cerebral hemodynamic movement 
which changes in association with an biological 
function, it is preferably displayed while being 
5 superimposed on the cerebral surface image. The 

cerebral surface image 23 is measured through three- 
dimensional MRI or three-dimensional X-ray CT and is 
displayed. The topographic image 22 is subjected to 
coordinate conversion such that coordinates of 

10 individual measuring points are positioned on the 
cerebral surface and values between individual 
measuring points subject to the coordinate conversion 
are interpolated to prepare the topographic image. 
When the prepared topographic image 22 is displayed 

15 while being superimposed on the cerebral surface image 
23, color of the overlying topographic image 22 is made 
to be semitransparent to allow the underlying cerebral 
surface image to be seen transparently. 

Figure 9 shows a diagram for explaining a 

20 measuring point coordinate converting method. A method 
of taking images of the form of three-dimensional MRI 
or three-dimensional X-ray CT will hereunder be 
described specifically. In Fig. 9, 5a designates an 
optical fiber for irradiation and 7a an optical fiber 

25 for detection. Information at a portion on the center 
line between these fibers is presumed for use as 
information at a desired measuring point. This is 
because a portion is used at which the supply of the 
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quantity of light from the irradiation fiber is 
maximized and the signal from an object to be measured 
is maximized. By performing imaging while disposing a 
marker at the measuring point 29 which is presumptively 
5 set by the living body optical measurement system, a 
skin and skeleton image 24, a cerebral image 25 and a 
marker image 26 can be displayed on the basis of the 
imaged form information. The images picked up as above 
have three-dimensional coordinate information. Thus, a 

10 perpendicular 28 passing through a measuring point 27 
indicated by the marker image 26 and being vertical to 
the skin surface at the measuring point 27 or the 
bottom of the marker image 26 is calculated and a point 
at which the perpendicular intersects the cerebral 

15 image 25 is defined as the measuring point 29 subject 
to the coordinate conversion. As shown in the present 
embodiment, it is known that when the cerebral function 
is measured, a change in hemodynamic movement having 
correlation to load occurs mainly at the cerebral 

20 surface (cerebral cortex). For the reasons set forth 
as above, by using anatomical information, the depth 
for coordinate conversion of measuring points can be 
known. But when an object to be measured is another 
living body organ such as muscles, the depth for 

25 coordinate conversion of measuring points cannot some- 
times be known from the anatomical information. In 
applying the present method to the measurement as 
above, light propagation of living tissue is calculated 



in advance through numerical calculation based on the 
Monte Carlo method to determine a depth best contri- 
buting to measured signals and coordinates of the 
measuring point are converted to the thus determined 
5 depth. 

The topographic image has been described as 
presumptive display image but a different image display 
method is available. 

For example, a square pixel of arbitrary size 

10 having the center of gravity at a presumptive measuring 
point is set at each presumptive measuring point, each 
pixel is displayed either in the form of an image in 
light and shade painted in color in correspondence to a 
value of each presumptive measuring point, the corre- 

15 spondence being determined in advance, or in the form 
of a bar graph image indicated by a bar or a length of 
line which corresponds to a value of each measuring 
point . 

In all image display methods, color arrange- 
20 ment can be selected freely when color is used but for 
measurement of hemodynamic movement, display in combi- 
nation of light and shade of red color and light and 
shade of blue color is preferable. This is because 
such an image that arterial hemodynamic is red and 
25 venous hemodynamic is blue is fixed. For example, the 
magnitude of positive measured value is displayed in 
light and shade of red color and the magnitude of 
negative measured value is displayed in light and shade 
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of blue color. 

(Embodiment 2) 

A second embodiment according to the present 
invention will be described. 
5 Fig. 10 shows an example of display of 

measured signals and estimated non-load signals . 
Displayed measured signals 110a and 110b are output 
signals from the lock-in amplifier 9a and estimated 
non-load signals 111a and 111b are calculated from the 

10 respective measured signals (calculation method will be 
described later) . 

The estimated non-load signals 111a and 111b 
are displayed on the display unit. In the displayed 
graph, abscissa represents measuring time and ordinate 

15 represents relative values of measured signals indi- 
cative of transmitting light intensity levels measured 
by the living body optical measurement system. 

When a load is applied to a subject, a 
loading start mark 112 indicative of a load application 

20 starting time point and a loading end mark 113 indi- 
cative of a load application ending time point are 
displayed in the form of straight lines. In the 
present embodiment, the cerebral cortex region dominat- 
ing the right-hand motion is measured from the scalp 

25 through the skull and the right-hand or left-hand 

motion is applied as a load (loads 1 and 3 correspond 
to the right-hand motion and loads 2 and 4 correspond 
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to the left-hand motion). Although in Fig. 10 all 
signals occurring throughout the measuring time, only a 
desired time interval (for example, a time interval 
including times before and after the loading time) can 
5 be displayed easily. Also, if the estimated non-load 
signals 111a and 111b are each displayed until a 
desired time point on an extension of the curve 
indicative of a time-variable change, the measured 
signals 110a and 110b and the estimated non-load 

10 signals 111a and 111b can be displayed simultaneously 
on real time base during measurement. By simultane- 
ously displaying the measured signals 110a and 110b and 
the estimated non-load signals 111a and 111b in this 
manner, a change in hemodynamic movement occurring in 

15 the living body can be decided easily by an observer. 
The estimated non-load signals displayed in advance on 
real time base can be corrected for display in the 
phase at which calculation of the estimated non-load 
signals is settled. 

20 The estimated non-load signals 111a and 111b 

can be determined by removing signals occurring during 
the load applying time (loading time) and signals 
occurring during the time for the signal to recover its 
intact form following removal of the load (relaxation 

25 time) and by fitting an arbitrary function to signals 
occurring during the remaining period through the 
method of least squares. Here, the arbitrary function 
and the relaxation time change with the kind of load 



and measurement locations and so those meeting the 
purpose of measurement are inputted through the input 
unit. In the present embodiment, the arbitrary 
function in the form of a polynomial of degree five is 
5 handled and the relaxation time is set to 30 seconds. 
Further, for the convenience of watching by the 
observer, different kinds of color or different kinds 
of lines can be used for display of signals. 

Figures 11A and 11B show examples of display 

10 of a difference signal between a measured signal and an 
estimated non-load signal and in these figures, a 
waveform of a difference signal 114a obtained by 
calculating the difference between the measured signal 
110a and the estimated non-load signal 111a in Fig. 10 

15 and a waveform of a difference signal 114b obtained by 
calculating the difference between the measured signal 
110b and the estimated non-load signal 111b in Fig. 10 
are displayed on the display unit. In the displayed 
graph, abscissa represents measuring time and ordinate 

20 represents relative difference signal intensity. 

Further, when a load is applied to a subject, a load 
start mark 112 indicative of a load application 
starting time point and a load end mark 113 indicative 
of a load application ending time point are displayed 

25 in the form of straight lines. Also, the present graph 
is a graph having its center at 0 and so shows a base 
line 115. 

In the present embodiment, waveforms 114a and 
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114b are displayed on different coordinate axes for 
different light source wavelengths but they can be 
displayed overlapping each other on the same coordinate 
axis. Also, for the convenience of watching by an 
5 observer, different kinds of color or different kinds 
of lines can be used for display. 

Figures 12A and 12B show examples of display 
of graphs depicting relative change amounts of concen- 
tration values of Hb02 and Hb (hereinafter represented 

10 by ACoxy and ACdeoxy, respectively) under the applica- 
tion of load. A waveform of a ACoxy signal 116a 
obtained from the measured signal 110a and estimated 
non-load signal 111a in Fig. 10 pursuant to equation 
(5) and a waveform of a ACdeoxy signal 116b obtained 

15 from the measured signal 110b and estimated non-load 
signal 111b in Fig. 10 pursuant to equation (5) are 
displayed on the display unit. In the displayed 
graph, abscissa represents measuring time and ordinate 
represents values of ACoxy and ACdeoxy. Further, a 

20 load start mark 112, a load end mark 113 and a base 
line 115 are also displayed. In the present embodi- 
ment, all intervals of measuring time are displayed but 
only desired time intervals (for example, a period 
including time points before and after the loading 

25 time) can be displayed. Here, the waveforms 116a and 
116b are displayed separately on different coordinate 
axes but they may be displayed overlapping each other 
on the same coordinate axis. Further, the individual 
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signals may be displayed in different kinds of color or 
in the form of different kinds of lines and for 
intuitive understanding by an observer, the ACoxy 
signal 116a may be displayed in, for example, a kind of 
5 red color and the ACdeoxy signal 116b may be displayed 
in, for example, a kind of green color. According to 
the measuring method and display method of the present 
invention, the correlation between load and measured 
signal is easy to understand and fluctuation is removed 

10 from the measured signal, thus ensuring that accuracy 
of signals can be increased. 

Figure 13 shows an example of display of a 
ACoxy loading time integral value 117a and a ACdeoxy 
loading time integral value 117b which are obtained 

15 during respective loading times. The ACoxy signal 114a 
and the ACdeoxy signal 114b in Figs. 11A and 11B are 
time integrated during each loading time to determine 
the ACoxy loading time integral value 117a and ACdeoxy 
loading time integral value 117b, which are displayed 

20 in the form of a cubic bar graph for respective load 
number. Here, abscissa represents load number and 
ordinate represents ACoxy loading time integral value 
and ACdeoxy loading time integral value. Alternative- 
ly, a ACoxy loading time average value and a ACdeoxy 

25 loading time average time may be displayed. Also, for 
the convenience of watching by an observer, display in 
different colors can be used. 

Figure 14 shows an example of display when 
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measurement is conducted at a plurality of measuring 
positions by using the living body optical measurement 
system. Here, an instance will be described where the 
portion to be measured is the head and four measuring 
5 positions are set on the head. 

In the present display example, a measuring 
portion image 118 of a subject, measuring position 
marks 119a to 119d representative of set measuring 
positions, graphs 121a to 121d corresponding to the 

10 respective measuring positions and index lines 120a to 
120d for indicating the correspondence relation between 
the respective measuring positions and the respective 
graphs are displayed on the display unit. Here, a head 
model figure or a measuring portion tomographic image 

15 or measuring portion three-dimensional image of a 

subject itself imaged by an image diagnostic apparatus 
represented by an MRI apparatus can be used as the 
measuring portion image 118. 

[Embodiment 3] 

20 The schematic construction of embodiment 3 of 

the living body optical measurement system according to 
the present invention is shown in Fig. 15. 

Light emitted from a light source 201 is 
collected using a lens system so as to impinge on an 

25 optical fiber 202 for light source. The light emitted 
from the light source is modulated in intensity with a 
desired frequency f of about 100Hz to 10MHz by means of 
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an oscillator 223 in order to remove noises due to 
external disturbance. Since the light source optical 
fiber 202 is connected to an optical fiber 204 for 
light irradiation through an optical fiber coupler 
5 203a, the light from the light source is transmitted to 
the light irradiation optical fiber 204 and is 
irradiated on a subject 206 by way of a light irradi- 
ation position 205. The wavelength of the light used 
depends on spectroscopic characteristics of an in vivo 

10 substance of interest but when the oxygen saturation 
amount and the blood amount are measured from concen- 
tration values of Hb and Hb02, a single or a plurality 
of wavelengths can be selected, for use, from light 
having the wavelength range of from 600nm to 1400nm. A 

15 semiconductor laser, a titanium/ sapphire laser or a 

light emitting diode can be used as the light source. 

Two light detection optical fibers 207a and 
^2 07b for detecting light transmitting through the 
subject 206 and going out of it are disposed at two 

20 different sites on the subject 206. In the present 

embodiment, the two light detection optical fibers 207a 
and 207b are disposed at two sites which are point 
symmetrical to a symmetry center of the light irradi- 
ation position 205. The light irradiation optical 

25 fiber 204 and the light detection optical fibers 207a 

and 207b are held in place by means of an optical fiber 
fixing member 208 having its surface painted in black. 
For simplicity of construction, the light irradiation 
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optical fiber 204, light detection optical fibers 207a 
and 207b and optical fiber fixing member 208 are 
integrally formed into a light detection probe to be 
detailed later. Since the light detection optical 
5 fibers 207a and 207b are connected to optical fibers 
209a and 209b for photodetectors through optical fiber 
couplers 203b and 203c, transmitting light rays 
detected by the light detection optical fibers 207a and 
207b are transmitted to photodetectors 210a and 210b 

10 and subjected to photoelectric conversion by the 

photodetectors 210a and 210b, so that transmitting 
light intensity levels are delivered in the form of 
electric signal intensity levels. Used as the photo- 
detectors 210a and 210b are photoelectric conversion 

15 elements such as for example photodiodes or photo- 
multiplier tubes . 

Of electric signals indicative of trans- 
mitting light intensity levels delivered out of the 
photodetectors 210a and 210b, only frequency components 

20 for light intensity modulation of the light source are 
extracted by lock-in amplifiers 224a and 224b, respec- 
tively. While an output from the lock-in amplifier 
224a is subjected to logarithmic conversion by a 
logarithmic amplifier 225a and then inputted to the 

25 negative pole of a differential amplifier 211, an 

output from the lock-in amplifier 224b is subjected to 
logarithmic conversion by a logarithmic amplifier 225b 
and subsequently inputted to the positive pole of the 
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differential amplifier 211. As a result, a difference 
signal between transmitting light intensity levels at 
the two different sites is delivered out of the 
differential amplifier 211 as an output signal. The 
5 output signal from the differential amplifier 211 is 
sequentially converted into a digital signal by an A/D 
converter 212, fetched into a computer 213 and 
displayed on a display unit 214 as time series data. 

Here, when a region 215 where the hemodynamic 

10 movement changes locally is included in only a view 
field 216b of the light detection optical fiber as 
shown in Fig. 15, the measured logarithmic difference 
signal reflects only a change in hemodynamic movement 
at the local region 215. On the presupposition that 

15 for near infrared rays, hemoglobin serving as a main 
constituent in hemodynamic dominantly acts on extinc- 
tion, the meaning of the measured logarithmic dif- 
ference signal will be described below. 

Where measuring time is t, light source 

20 wavelength is X, irradiation light intensity is 10 (t), 
concentration values of hemoglobin oxide and reduced 
hemoglobin are Cox(t) and Cdeox(t), respectively, 
changes in concentration values of hemoglobin oxide and 
reduced hemoglobin which occur at the local region 215 

25 are ACox(t) and ACdeox(t), respectively, extinction 
coefficients for light source wavelength X of hemo- 
globin oxide and reduced hemoglobin are €ox(A) and 
Gdeox(JL), respectively, attenuation due to scattering 
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and absorption caused by other constituents than 
hemoglobin is Ds and weight coefficient caused by 
scattering is d, transmitting light intensity signal 
Id(t) detected by the photodetector 210b is given by 
the following equation (6) and transmitting light 
intensity signal Id'(t) detected by the photodetector 
210a is given by the following equation (7): 

Id(t)=Ds. exp[-[€ox(X) (Cox( t)+ACox(t ) ) 

+€deox(A.) <Cdeox(t)+ACdeox(t) ) ]I0(t) (6) 
Id' (t)=Ds. exp[-[eox(A)Cox(t)+Gdeox(A.)Cdeox(t) ]d]l0(t) 

(7) 

Next, equations (6) and (7) are expressed in 
terras of natural logarithm and then equation ( 7 ) is 
subtracted from equation ( 6 ) to obtain the following 
equation ( 8 ) . The left side of equation ( 8 ) is the 
measured logarithmic difference signal. 

ln[Id(t)/Id' (t)]=-[Gox(X)ACox(t)+€deox(A)ACdeox(t)]d 

(8) 

Here, when 805nm±10nm is particularly used as 
the light source wavelength for measurement, 

Gox(805±10)^€deox( 805+10) (9) 



stands and by using constant K, equation (8) can be 
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reduced to 

ln[Id(t)/Id' (t) ]=-[ACox(t)+ACdeox(t) ]. K 

(10) 

Accordingly, the logarithmic difference signal 
5 measured by using the light source wavelength 805nm± 
lOnm represents a value corresponding to a change 
amount in hemodynamic amount [ ACox ( t ) + ACdeox ( t ) ] 
(hereinafter referred to as relative hemodynamic change 
amount) . Also, when the number of wavelengths used for 

10 the light source is set to two (,Xl,A.2), different 

intensity modulation frequencies (fl,f2) are applied to 
the respective wavelengths and frequency separation is 
effected by means of the lock-in amplifiers, trans- 
mitting light intensity signals of the individual 

15 wavelengths can be measured. Accordingly, equation (8) 
holds for the respective wavelengths and a simultaneous 
equation consisting of the following equations (11) and 
(12) can be introduced: 

ln[Id(Xl, t)/Id'(A.l, t)] 
20 =-[6ox(Al)ACox(t)+€deox(Al)ACdeox(t) ]d (11) 

ln[Id(A2, t)/Id'(A.2, t) 
=-[Gox(X2)ACox(t)+Gdeox(X2)ACdeox(t) ]d (12) 



Since the extinction coefficients eox(Al), 
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€ox(A2), edeox(Al) and 6deox(A2) are known, value 
ACox(t)d corresponding to a change amount of hemoglobin 
oxide and value ACdeox(t) corresponding to a change 
amount of reduced hemoglobin can be determined by 
5 solving equations (11) and (12) with the computer 213 
and time series data representative of a determined 
relative change amount can be displayed graphically on 
the display unit 214. For expansion of the above, the 
number of wavelengths can be increased, d can be erased 

10 or a relative change amount of concentration of a light 
absorbing substance, other than hemoglobin, which 
exists in a small amount can be determined. 

Alternatively, as shown in Fig. 16, the lock- 
in amplifiers, logarithmic amplifiers and differential 

15 amplifiers may not be used, detection signals from the 
photodetectors 210a and 210b may be converted into 
digital signals by A/D converters 212, respectively, 
the digital signals may then be FFT processed with the 
computer 213 to extract only a signal corresponding to 

20 the intensity modulation frequency of the light source, 
a logarithmic difference between transmitting light 
intensity levels at two different detection positions 
is calculated through a similar procedure to the above 
calculation process, and the determined relative change 

25 amount can be displayed graphically as time series data 
on the display unit 214. 

Figures 17A and 17B show an example of the 
light detection probe. Figure 17A illustrates a 
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section of the light detection probe and Fig. 17B is a 
diagram of the light detection probe as viewed from a 
subject contact surface. 

The light detection probe is comprised of the 
5 single light irradiation optical fiber 204, the two 
light detection optical fibers 207a and 207b and the 
optical fiber fixing member 208 having its surface 
painted in black and made of metal or plastics, and the 
optical fibers are connected with the optical fiber 

10 couplers 203a, 203b and 203c. In order to maintain 
flexibility of the individual optical fibers, each 
optical fiber is constructed of a plurality of optical 
fibers. As a material of the optical fiber, plastics 
or quartz is used. When the present light detection 

15 probe is used for a living body, the subject contact 
surface 217 is covered with, for example, resilient 
sponge . 

The size of the detection surface of each 
light detection optical fibers 207a or 207b must be 

20 changed in accordance with the purpose and the state of 
a subject but when measurement of, for example, the 
cerebral function is performed, the sectional shape is 
made to be a circle having a diameter of about 1mm to 
20mm or a square having a side of about 1mm to 20mm. 

25 The two light detection optical fibers 207a and 207b 
are disposed at positions r (r=5mm to 50mm) distant 
from the light irradiation optical fiber 204 and here, 
disposed symmetrically. A plurality of kinds of light 
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detection probes, in which the distance r differs and 
the sectional shape of each light detection optical 
fiber 207a or 207b differs, are prepared and one probe 
is exchanged with another in compliance with the 
5 purpose of measurement, thereby permitting convenient 
measurement. Since the light reaching depth substan- 
tially equals the distance r from the light source, a 
depth approximating the cerebral cortex of the cerebrum 
can be measured from the head surface through the 
10 skull. 

In the light detection probe, various forms 
of disposition of the light detection optical fibers 
207 can be conceivable. For example, as shown in Fig. 
18, four light detection optical fibers 207a, 207b, 

15 207c and 207d can be disposed at positions r 

equidistant from the light irradiation optical fiber 
204, and desired two light detection optical fibers can 
be selected for measurement. In an alternative, any 
optical fiber may not be used but a lens system may be 

20 used or a light source and photodetectors can be 
disposed directly in the fixing member 208. 

Figure 19 shows an example where the optical 
measuring system according to the present invention is 
used for measurement of the cerebrum of a living body. 

25 A light detection probe comprised of optical fiber 
couplers 203a, 203b and 203c, a light irradiation 
optical fiber 204, light detection optical fibers 207a 
and 207b, and an optical fiber fixing member 208 is 
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fixed to a subject 206 by means of a fitting belt 218 
made of rubber. The light irradiation optical fiber 
204 is connected to a light source optical fiber 202 
through the optical fiber coupler 203a and the light 
5 detection optical fibers 207a and 207b are connected to 
optical fibers 209a and 209b for light detection, 
respectively, through the optical fiber couplers 203b 
and 203c. Provided on a front panel of the optical 
measuring system 219 are connectors for the light 

10 source optical fiber 202 and light detection optical 
fibers 209a and 209b, an output signal adjusting knob 
220, an output signal value indicating window 221 and a 
display unit 214. Arranged in the optical measuring 
system 219 are the differential amplifiers, A/D 

15 converters, microprocessor, light source, photo- 
detectors, optical switches and other necessary 
electric circuits . 

A logarithmic difference signal value between 
transmitting light intensity levels detected at two 

20 sites is digitally indicated on the output signal value 
indicating window 221 and an offset value of the 
logarithmic difference signal value is determined using 
the output signal adjusting knob 220. For example, in 
the absence of a local change in hemodynamic movement 

25 in the cerebrum of the subject, the logarithmic dif- 
ference signal between transmitting light intensity 
levels detected at the two sites is so adjusted as to 
be zero. Thereafter, measurement is started and time 
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series data 222 representative of the logarithmic 
difference signal is graphically displayed on the 
display unit 214. Also, the arithmetic operation 
described previously is carried out to graphically 
5 display a local hemodynamic amount or time-variable 

changes in relative change amounts of hemoglobin oxide 
quantity and reduced hemoglobin quantity. 

[Embodiment 4] 

The schematic construction of embodiment 4 of 

10 the living body optical measurement system according to 
the present invention is shown in Fig. 20. 

Light emitted from a light source 201 is 
collected using a lens system so as to impinge on an 
optical fiber 202 for light source. The light emitted 

15 from the light source is modulated in intensity with a 
desired frequency of about 100Hz to 10MHz by means of 
an oscillator 223 in order to remove noises due to 
external disturbance. Since the light source optical 
fiber 202 is connected to an optical fiber 204 for 

20 light irradiation through an optical fiber coupler 

203a, the light from the light source is transmitted to 
the light irradiation optical fiber 204 and is irradi- 
ated on a subject 206 by way of a light irradiation 
position 205. The wavelength of the light used depends 

25 on spectroscopic characteristics of an in vivo sub- 
stance of interest but when the oxygen saturation 
amount and the blood amount are measured from 
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concentration values of Hb and Hb02, a single or a 
plurality of wavelengths can be selected, for use, from 
light having the wavelength range of from 600nm to 
1400nm. A semiconductor laser, a titanium/sapphire 
5 laser or a light emitting diode can be used as the 
light source . 

Four light detection optical fibers 207a, 
207b, 207c and 207d for detecting light transmitting 
through the subject 206 and going out of it are 

10 disposed at four different sites on the subject 206. In 
the present embodiment, the two light detection optical 
fibers 207b and 207c are disposed at two sites which 
are point symmetrical to a symmetry center of the light 
irradiation position 205, the light detection optical 

15 fiber 207a is disposed such that the centroid point of 
the light detection optical fiber 207a is on a half- 
line having its origin at the centroid point of the 
light irradiation position and passing through the 
centroid point of the light detection optical fiber 

20 207b, and the light detection optical fiber 207d is 
disposed such that the centroid point of the light 
detection optical fiber 207d is on a half-line having 
its origin at the centroid point of the light irradi- 
ation position and passing through the centroid point 

25 of the light detection optical fiber 207c. The light 
detection optical fibers 207a and 207d can be disposed 
anywhere so long as their centroid points are on the 
half -lines but in the present embodiment, they are 
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disposed point-symmetrically to the symmetry center of 
the light irradiation position 205 and outside the 
light detection optical fibers 207b and 207c. Here, 
the light irradiation optical fiber 204 and the light 
5 detection optical fibers 207a, 207b, 207c and 207d are 
held in place by means of an optical fiber fixing 
member 208 made of metal and having its surface painted 
in black. Since the light detection optical fibers 
207a, 207b, 207c and 207d are connected to optical 

10 fibers 209a, 209b, 209c and 209d for photodetectors 
through optical fiber couplers 203b, 203c, 203d and 
203e, transmitting light rays detected by the light 
detection optical fibers 207a, 207b, 207c and 207d are 
transmitted to photodetectors 210a, 210b, 210c and 210d 

15 and subjected to photoelectric conversion by the photo- 
detectors 210, so that transmitting light intensity 
levels subject to the photoelectric conversion are 
delivered in the form of electric signal intensity 
levels. Used as the photodetectors 210 are photo- 

20 electric conversion elements such as for example 
photodiodes or photomultiplier tubes. 

Of electric signals indicative of trans- 
mitting light intensity levels delivered out of the 
photodetectors 210a and 210b, only frequency components 

25 for intensity modulation of the light source are 

extracted by lock— in amplifiers 224a and 224b, respec- 
tively. While an output from the lock-in amplifier 
224a is subjected to logarithmic conversion by a 
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logarithmic amplifier 225a and then inputted to the 
negative pole of a differential amplifier 211a, an 
output from the lock- in amplifier 224b is subjected to 
logarithmic conversion by a logarithmic amplifier 225b 
5 and subsequently inputted to the positive pole of the 
differential amplifier 211a. Of electric signals 
indicative of transmitting light intensity levels 
delivered out of the photodetectors 210c and 210d, only 
frequency components for intensity modulation of the 

10 light source are extracted by lock-in amplifiers 224c 

and 22 Ad, respectively. While an output from the lock- 
in amplifier 224d is subjected to logarithmic conver- 
sion by a logarithmic amplifier 225d and then inputted 
to the negative pole of a differential amplifier 211b, 

15 an output from the lock-in amplifier 224c is subjected 
to logarithmic conversion by a logarithmic amplifier 
225c and subsequently inputted to the positive pole of 
the differential amplifier 211b. Further, an output 
from the differential amplifier 211a is inputted to the 

20 negative pole of a differential amplifier 211c and an 

output from the differential amplifier 211b is inputted 
to the positive pole of the differential amplifier 
211c. As a result, a difference signal among trans- 
mitting light intensity levels at the four different 

25 sites is delivered out of the differential amplifier 
211c as an output signal. The output signal from the 
differential amplifier 211c is sequentially converted 
into a digital signal by an A/D converter 212, fetched 
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into a computer 213 and displayed graphically on a 
display unit 214 as time series data. 

Here, when a region 215 where the hemodynamic 
movement changes locally is included in only a view 
field 216b of the light detection optical fiber as 
shown in Fig. 20, the logarithmic differential signal 
of transmitting light intensity delivered out of the 
differential amplifier 211c reflects only a change in 
local hemodynamic movement. On the presupposition that 
for near infrared rays, hemoglobin serving as a main 
constituent in hemodynamic dominantly acts on extinc- 
tion, the meaning of the logarithmic difference signal 
delivered out of the differential amplifier 211c will 
be described below. 

Where measuring time is t, light source 
wavelength is X, irradiation light intensity is I0(t), 
concentration values of hemoglobin oxide and reduced 
hemoglobin are Cox(t) and Cdeox(t), respectively, 
changes in concentration values of hemoglobin oxide and 
reduced hemoglobin which occur at the local region 215 
are ACox(t) and ACdeox(t), respectively, extinction 
coefficients for light source wavelength k of hemo- 
globin oxide and reduced hemoglobin are 6ox(X) and 
Gdeox(X), respectively, attenuation contained in 
transmitting light intensity levels detected by the 
photodetectors 210b and 210c and due to scattering and 
absorption caused by other constituents than hemoglobin 
is Dsl, attenuation contained in transmitting light 
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intensity levels detected by the photodetectors 210a 
and 210d and due to scattering and absorption caused by 
other constituents than hemoglobin is Ds2, weight 
coefficient contained in the transmitting light 
5 intensity levels detected by the photodetectors 210b 
and 210c and caused by scattering is dl, and weight 
coefficient contained in the transmitting light 
intensity levels detected by the photodetectors 210a 
and 210d and caused by scattering is d2, transmitting 

10 light intensity signal Idl(t) detected by the photo- 
detector 210c r transmitting light intensity signal 
Id2(t) detected by the photodetector 21 Od, transmitting 
light intensity signal Idl'(t) detected by the photo- 
detector 210b and transmitting light intensity signal 

15 Id2'(t) detected by the photodetector 210a are given by 
the following equations (13) to (16): 

Idl(t)=Dsl. exp[-[€ox(A) (Cox(t)+ACox(t) ) 

+edeox(X) (Cdeox(t)+ACdeox(t) ) ]dl]l0(t) 

(13) 

20 Id2(t)=Ds2. expt-t€ox(A.)(Cox(t)+ACox(t) ) 

+edeox(A.) (Cdeox(t)+ACdeox(t) ) ]d2]I0(t) 

(14) 

Idl' (t)=Dsl. exp[-[Gox(A.)Cox(t) 

+€deox(A)Cdeox(t) ]dl]I0(t) (15) 

25 Id2'(t)=Ds2. exp[-[eox(A)Cox(t) 

+edeox(A)Cdeox(t)]d2]I0(t) (16) 
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Next, equations (13) and (14) are expressed 
in terms of natural logarithm and then equation (14) is 
subtracted from equation (13) to obtain the following 
equation (17): 

ln[ Idl ( t ) /Id2 ( t ) ] =ln[ Dsl/Ds2 ] 
-[6ox(X) (Cox(t)+ACox(t) ) 
+€deox ( X ) ( Cdeox ( t ) +ACdeox ( t ) ) ] ( dl-d2 ) 

(17) 

Equations (15) and (16) are expressed in 
terms of natural logarithm and then equation (16) is 
subtracted from equation (15) to obtain the following 
equation (18): 

ln[Idl' (t)/Id2' (t) ]=ln[Dsl/Ds2] 
-Eox(X) (Cox(t)+ACox(t) ) 
+6deox ( X ) ( Cdeox ( t ) +ACdeox ( t ) ) ] ( dl-d2 ) 

(18) 

The left side of equation (17) represents the 
output of the differential amplifier 211b and the left 
side of equation (18) represents the output of the 
differential amplifier 211a. Here, by subtracting 
equation (18) from equation (17), the following 
equation (19) is obtained: 
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ln[ (Idl(t)/Id2(t) )(Id2' (t)/ldl' (t)] 

=-[€ox(A)ACox(t)+€deox(A)ACdeox(t) ] (dl-d2) 

(19) 

The left side of equation (19) represents the 
output of the differential amplifier 211c, that is, the 
measured logarithmic difference signal. 

Here, when 805nm±10nm is particularly used as 
the light source wavelength for measurement, the afore- 
mentioned relation of equation ( 9 ) stands and by using 
constant K, equation (19) can be reduced to the 
following equation (20): 

ln[ (Idl(t)/Id2(t) ) (Id2' (t)/Idl' (t) ) ] 

=-[ACox(t)+ACdeox(t) ]. K (20) 

Accordingly, the logarithmic difference 
signal measured using the light source wavelength 
805nm±10nm represents a value corresponding to a 
relative hemodynamic change amount [ ACox ( t ) +ACdeox ( t ) J . 

Also, when the number of wavelengths used for 
the light source is two (Xl,k2) r different intensity 
modulation frequencies (fl,f2) are applied to the 
respective wavelengths and frequency separation is 
effected by means of the lock-in amplifiers, trans- 
mitting light intensity signals of the individual 
wavelengths can be measured. Accordingly, equation 
(19) holds for the respective wavelengths and a 
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simultaneous equation consisting of the following 
equations (21) and (22) can be introduced: 

ln[(Idl(Al, t)/Id2(Al, t))(Id2'(Al, t)/Idl'(Al, t))] 
=-[6ox(Xl)ACox(t)+Gdeox(Xl)ACdeox(t) ] (dl-d2) 
5 (21) 

ln[Idl(A2, t)/Id2(A2, t))(Id2'(A2, t)/Idl'(A2, t))] 
=-[€ox(A2)ACox(t)+Gdeox(A2)ACdeox(t) ] (dl-d2) 

(22) 

Since the extinction coefficients €ox(Al), 
10 €ox(A2), edeox(Al) and €deox(A2) are known, value 
ACox( t) (dl-d2 ) corresponding to a change amount of 
hemoglobin oxide and value ACdeox(t) (dl-d2) correspond- 
ing to a change amount of reduced hemoglobin can be 
determined by solving equations (21) and (22) with the 
15 computer 213 and time series data representative of a 
determined relative change amount can be displayed 
graphically on the display unit 214. For expansion of 
the above, the number of wavelengths can be increased, 
(dl-d2) can be erased or a relative change amount of 
20 concentration of a light absorbing substance, other 

than hemoglobin, which exists in a small amount can be 
determined . 

Alternatively, as shown in Fig. 21, the lock- 
in amplifiers, logarithmic amplifiers and differential 
25 amplifiers may not be used, detection signals from the 
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photodetectors 210a, 2 10b f 210c and 210d may be 
converted into digital signals by the A/D converters 
212, respectively, the digital signals may then be FFT 
processed with the computer 213 to extract only a 
5 signal corresponding to the intensity modulation 

frequency of the light source, a logarithmic difference 
among transmitting light intensity levels at four 
different detection positions may be calculated through 
a similar procedure to the above calculation process, 
10 and then the determined relative change amount can be 
displayed graphically as time series data on the 
display unit 214. 

[Embodiment 5] 

The schematic construction of embodiment 5 of 

15 the living body optical measurement system according to 
the present invention is shown in Fig. 22. 

Light rays emitted from light sources 201a 
and 201b are collected using lens systems so as to 
impinge on an optical fiber 202a for light source and 

20 an optical fiber 202b for light source, respectively. 
The light rays emitted from the respective light 
sources are modulated in intensity with desired dif- 
ferent frequencies f of about 100Hz to 10MHz by means 
of oscillators 223a and 223b in order to remove noises 

25 due to external disturbance. Here, the intensity 

modulation frequency for the light source 201a is set 
to f 1 and the intensity modulation frequency for the 
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light source 201b is set to f2. Since the light source 
optical fiber 202a is connected to an optical fiber 
204a for light irradiation through an optical fiber 
coupler 203a and the light source optical fiber 202b is 
5 connected to an optical fiber 204b for light irradi- 
ation through an optical fiber coupler 203c, the light 
rays from the respective light sources are transmitted 
to the light irradiation optical fibers 204a and 204b 
and irradiated on a subject 206 by way of light irradi- 

10 ation positions 205a and 205b. In order to obtain 

reference light rays, the light is split at a halfway 
point of each of the light irradiation optical fibers 
204a and 204b by means of a splitter 226a or 226b and 
the intensity of each light source is converted into an 

15 electric signal by means of a photodetector 210a or 

210c. A reference light intensity signal for the light 
source 201a, delivered out of the photodetector 210a, 
is inputted to a lock-in amplifier 224a and separated 
on the basis of a reference frequency from the oscil- 

20 lator 223a. An output of the lock-in amplifier 224a is 
inputted to a logarithmic amplifier 225a so as to 
undergo logarithmic conversion and then inputted to the 
negative pole of a differential amplifier 211a. A 
reference light intensity signal for the light source 

25 201b, delivered out of the photodetector 210c, is 

inputted to a lock-in amplifier 224d and separated on 
the basis of a reference frequency from the oscillator 
223b. An output of the lock- in amplifier 224d is 
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inputted to a logarithmic amplifier 225d so as to 
undergo logarithmic conversion and then inputted to the 
negative pole of a differential amplifier 211b. The 
wavelength of the light used depends on spectroscopic 
5 characteristics of an in vivo substance of interest but 
when the oxygen saturation amount and the hemodynamic 
amount are measured from concentration values of Hb and 
Hb02, a single or a plurality of wavelengths can be 
selected, for use, from light having the wavelength 

10 range of from 600nm to 1400nm. A semiconductor laser, 
a titanium/sapphire laser or a light emitting diode can 
be used as the light source. 

For detection of light transmitting through a 
subject 206 and going out of it, one optical fiber 207 

15 for light detection is disposed at a position on the 

subject 206 which is equidistant from light irradiation 
positions 205a and 205b. Here, the light irradiation 
optical fibers 204a and 204b and the light detection 
optical fiber 207 are held in place by means of an 

20 optical fiber fixing member 208 having its surface 

painted in black. Since the light detection optical 
fiber 207 is connected to an optical fiber 209 for 
light detection through an optical fiber coupler 203b, 
transmitting light detected by the light detection 

25 optical fiber 207 is transmitted to a photodetector 

210b and subjected to photoelectric conversion by the 
photodetector 210b, so that the transmitting light 
intensity is delivered in the form of electrical signal 



- 62 - 

intensity. Used as the photodetector 210b is a photo- 
electric conversion element such as for example a 
photodiode or a photomultiplier tube. 

Since the electric signal representative of 
5 the transmitting light intensity level delivered out of 
the photodetector 210b contains the transmitting light 
intensity signal for the light source 201a and the 
transmitting light intensity signal for the light 
source 201b, only an intensity modulation frequency 

10 component for the light source 201a is extracted by a 

lock-in amplifier 224b and only an intensity modulation 
frequency component signal for the light source 201b is 
extracted by a lock-in amplifier 224c. An output from 
the lock-in amplifier 224b is subjected to logarithmic 

15 conversion by a logarithmic amplifier 225b and then 
inputted to the positive pole of the differential 
amplifier 211a. An output from the lock-in amplifier 
224c is subjected to logarithmic conversion by a 
logarithmic amplifier 225c and subsequently inputted to 

20 the positive pole of the differential amplifier 211b. 
As a result, a logarithmic difference signal between 
the intensity of the light source 201a and the trans- 
mitting light intensity for the light source 201a is 
delivered out of the differential amplifier 211a as an 

25 output signal, and a logarithmic difference signal 

between the intensity of the light source 201b and the 
transmitting light intensity for the light source 201b 
is delivered out of the differential amplifier 211b as 



an output signal. Further, an output from the dif- 
ferential amplifier 211a is inputted to the negative 
pole of a differential amplifier 211c and an output 
from the differential amplifier 211b is inputted to the 
positive pole of the differential amplifier 211c, so 
that a logarithmic differential signal between trans- 
mitting light intensity levels which is removed of 
fluctuation in the light source intensity can be 
delivered out of the differential amplifier 211c. The 
output signal from the differential amplifier 211c is 
sequentially converted into a digital signal by an A/D 
converter 212, fetched into a computer 213 and dis- 
played on a display unit 214 as time series data. 

Here, when a region 215 where the hemodynamic 
movement changes locally is included in only a view 
field 216b of the light detection optical fiber as 
shown in Fig. 22, the measured logarithmic differential 
signal reflects only a change in local hemodynamic 
movement. On the presupposition that for near infrared 
rays, hemoglobin serving as a main constituent in 
hemodynamic dominantly acts on extinction, the meaning 
of the measured logarithmic difference signal will be 
described below. 

Where measuring time is t, light source 
wavelength is k, irradiation light intensity from the 
irradiation position 205b is 10 (t), irradiation light 
intensity from the irradiation position 205a is 10' (t), 
reference light intensity from the splitter 226b is 
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Ir(t), reference light intensity from the splitter 226b 
is Ir'(t), ratio of splitting to reference light of the 
splitter is <x indicating that 



5 I0(t):lr(t)=10'(t):lr'(t)=l:a, 

concentration values of hemoglobin oxide and reduced 
hemoglobin are Cox(t) and Cdeox(t), respectively, 
changes in concentration values of hemoglobin oxide and 
reduced hemoglobin which occur at the local region 215 

10 are ACox(t) and ACdeox, respectively, extinction 

coefficients for light source wavelength k of hemo- 
globin oxide and reduced hemoglobin are Gox(X) and 
edeox(A.), respectively, attenuation due to scattering 
and absorption caused by other constituents than 

15 hemoglobin is Ds, and weight coefficient caused by 

scattering is d, transmitting light intensity signal 
Id(t) for the light source 201b detected by the 
photodetector 210b, that is, the output from the lock- 
in amplifier 224c is given by the following equation 

20 (23) and transmitting light intensity signal Id'(t) for 
the light source 201a, that is, the output from the 
lock-in amplifier 224b is given by the following 
equation (24): 



25 



Id(t)=Ds. exp[-[Gox(A) (Cox(t)+ACox(t) 

+edeox(A) (Cdeox(t)+ACdeox(t) )]d]I0(t) 

(23) 



Id'(t) 
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=Ds. exp[-[Gox(A.)Cox(t) + 
Gdeox(A)Cdeox(t) ]d]IO' (t) 



(24) 



Next, equations (23) and (24) are expressed 
in terms of natural logarithm and then equation (23) is 
reduced to the following equation (25) and equation 
(24) is reduced to the following equation (26): 

ln[Id(t)/lO(t)]=ln[Ds]-[6ox(X) (Cox(t)+ACox(t ) 

+6deox ( A) ( Cdeox ( t ) +ACdeox ( t ) ) ] d 

(25) 

ln[Id' (t)/lO' (t)]=ln[Ds]-[Gox(A)Cox(t) 

+6deox(A)Cdeox(t) ]d (26) 

Further, equation (26) is subtracted from 
equation (25) to obtain the following equation (27): 

ln[(Id(t)/Id' (t)) (10' (t)/lO(t)] 
=-[€ox(A)ACox(t)+edeox(A)ACdeox(t)]d (27) 

Here, 

Ir(t)=al0(t) (28) 

Ir' (t)=ctlO' (t) (29) 

stand and therefore, the output from the differential 
amplifier 211a is 



ln[Id' (t)/alO' (t) ] 
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and the output from the differential amplifier 211c is 

ln[ <Id(t)/Id' (t) ) (10' (t)/I0(t)] (30) 

Since equation (30) equals the left side of equation 
5 (27), the logarithmic difference signal delivered out 
of the differential amplifier 211c is equivalent to 
equation (27). 

Here, when 805nm±10nm is particularly used as 
the light source wavelength for measurement, the afore- 
10 mentioned relation of equation (9) stands and by using 
constant K, equation (27) can be reduced to the 
following equation (31): 

In[(ld(t)/Id'(t))(l0'<t)/I0(t))] 
=-[ACox(t)+ACdeox(t) ]K (31) 

15 

Accordingly, the logarithmic difference 
signal measured using the light source wavelength 
805nm±10nm represents a value corresponding to a 
relative hemodynamic change amount [ ACox ( t ) +ACdeox ( t ) ] - 
20 K. Also, when the number of wavelengths used for the 
light source is two (Al,A.2), intensity modulation 
frequencies (fl, f2, f3, f4) which differ with the 
respective wavelengths and the respective irradiation 
positions are applied and frequency separation is 
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effected by means of the lock-in amplifiers, trans- 
mitting light intensity signals for the individual 
wavelengths and the respective irradiation positions 
can be measured. Accordingly, equation (27) holds for 
5 the respective wavelengths and a simultaneous equation 
consisting of the following equations (32) and (33) can 
be introduced: 



ln[ (Id(A.l,t)/Id' (Al,t) ) (10' (Al,t)/I0(Al,t) ) ] 
=- [ eox ( XI ) ACox ( t )+€deox ( A.1 ) ACdeox( t ) ]d 
10 (32) 

ln[ (Id(A.2,t)/Id' (A.2,t) ) (10' (A.2,t)/I0(X2,t) ) ] 
=- [ 6ox ( X2 ) ACox ( t ) +Gdeox ( X2 ) ACdeox ( t ) ]d 

(33) 



Since the extinction coefficients Gox(Al), 
15 6ox(X2), Gdeox(Al) and €deox(A2) are known, value 

ACox(t)d corresponding to a change amount of hemoglobin 
oxide and ACdeox(t) corresponding to a change amount of 
reduced hemoglobin can be determined by solving equa- 
tions (32) and (33) with the computer 213 and time 
20 series data representative of a determined relative 
change amount can be displayed graphically on the 
display unit 214. For expansion of the above, the 
number of wavelengths can be increased, d can be erased 
or a relative change amount of concentration of a light 
25 absorbing substance, other than hemoglobin, which 
exists in a small amount can be determined. 
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Alternatively, as shown in Fig. 23, the lock- 
in amplifiers, logarithmic amplifiers and differential 
amplifiers may not be used, detection signals from the 
photodetectors 210a, 210b and 210c may be converted 
into digital signals by the A/D converters 212, respec- 
tively, the digital signals may then be FFT processed 
with the computer 213 to extract only signals corre- 
sponding to intensity modulation frequencies of the 
respective light sources, a relative change amount 
calculated and determined through a similar procedure 
to the above calculation process can be displayed 
graphically as time series data on the display unit 
214. 

In the present invention, since the low cost 
light irradiation means and photodetectors are used to 
realize the simplified arithmetic processing, the high- 
speed processing can be ensured with an economical 
system and the living body function in correspondence 
to a plane image illustrative of the shape of an object 
to be measured can be imaged, thus ensuring that means 
effective especially for measurement of the local 
function of living body can be provided. 

According to the prior arts, when the 
hemodynamic movement of a subject is measured by 
alternately repeating unloading and loading, the 
reference value varies and correct measurement cannot 
be carried out unless the measurement is delayed until 
signals to be measured are stabilized by keeping the 
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subject quiet. According to the measuring method of 
the present invention, measurement can be conducted 
without waiting for the stabilization of signals . 
Further, measured signals can be removed of fluctuation 
5 to promote accuracy of signals . 

When the first and second detection positions 
are set at locations which are substantially equi- 
distant from the light irradiation position, trans- 
mitting light intensity signals at the respective 

10 detection positions equally change with an overall 
change in hemodynamic movement in a living body. 
Accordingly, by taking a logarithmic difference between 
a transmitting light intensity signal at the first 
detection position and a transmitting light intensity 

15 signal at the second detection position, a signal 

change attributable to the overall hemodynamic movement 
change can be removed. Further, when a change attri- 
butable to local hemodynamic movement is included in 
only one of transmitting light intensity signals at the 

20 first and second detection positions, the transmitting 
light intensity logarithmic difference signal reflects 
only the local change in hemodynamic movement. 

The light ray comes into the living body 
through the light irradiation position and passes 

25 through complicated paths while interacting with 
various tissues of the living body and undergoing 
scattering and attenuation before it goes out of the 
living body through the light detection position. In 
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the present invention, since the logarithmic difference 
between light intensity levels going out of the living 
body through positions which are substantially 
equidistant from the light irradiation position, the 
influence of scattering and attenuation caused by the 
tissues of the living body can be cancelled out and a 
slight signal reflecting the local change in hemo- 
dynamic movement can be detected with high accuracy. 

Concretely, when an operation for cutting off 
a portion of the cerebrum is conducted in determination 
of the focus of epilepsy and therapy of serious 
epilepsy, an image of the portion to be cut off can be 
used to confirm that the cutting-off portion is a 
portion which does not impair the important in vivo 
function of the subject. 



